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ABSTRACT 

Recent developments in pseudo-holography provide the means to 
generate a three-dimensional image of an object from any radiation it 
scatters. One promising technique uses scattered radiation to cast 
shadows of a Fresnel zone plate, referred to as the coded image, 
which are recorded on film. The coded image is then decoded by 
passing coherent light through it which is diffracted producing an 
image of the radiation source. 

In this paper, a brief illustration of pseudo-holography using 
scattered thermal neutrons is presented, along with a development of 
a general expression for the coded image. The remainder of the paper 
is devoted to developing a scattered neutron flux source function by 
diffusion theory methods for a non-reentrant, steady state, homogeneous 
body of simple geometry. To test the accuracy of the source function, 
the computed scattered flux is compared with some measured values. 

The comparison supports the diffusion approximation for detection 
methods which do not have a strong forward detecting preference. 
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I. INTRODUCTION 

Recent developments in pseudo-holography provide the means to 

generate a three-dimensional image of an object from any radiation it 

scatters. One promising pseudo-holographic technique is basically a 

two-step process. The scattered radiation is used to cast shadows of 

a Fresnel zone plate (hereafter termed the coded image) which are 

recorded on film. The coded image is then decoded by passing coherent 

light through it which is diffracted producing an image of the radiation 

scattering source. This technique, first proposed for radiation sources 

by Mertz and Young^ in x-ray astronomy, and extended to gamma ray 

2 3 A’ 5 6 

sources by Barrett * * * and Rogers , has been successfully applied to 

actual image formation using x-rays scattered from an object by Sider, 

et al.^ The use of scattered neutrons in image formation is presently 

0 

under study by Rose, et al. and provides the motivation for this paper. 

The pseudo-holographic technique using scattered neutrons differs 
from the use of scattered x-rays only in that the film and Fresnel zone 
plate respond to neutrons, not x-rays, in forming the coded image. It 
follows that an analytical representation of the coded image can be 
developed by describing the scattered neutron flux, modifying the 
scattered flux with a Fresnel zone plate response, and applying a film 
(detector) response function to the modified flux (shadows) seen by the 
film. The first step in this analysis is an accurate description of the 
uninterrupted scattered flux. It is the intent of this paper to develop 
an analytical method which describes scattered neutron fluxes emerging 
from objects in a beam of thermal neutrons. 
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First, a brief illustration of pseudo-holography is presented 
along with a development of a general expression for the coded image. 
The remainder of the paper is devoted to developing a scattered flux 
source function by diffusion theory methods for a non-reentrant, 
steady state, homogeneous body of simple geometry. To test the 
accuracy of the source function the computed scattered flux is compared 
with some measured values. 
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II. BACKGROUND 



A« Concept Illustration 

The basic concept of the pseudo-holographic technique noted above, 
as applied to thermal neutrons, is illustrated in Figure 1. A Fresnel 
zone plate is positioned relative to the scattering body and a suitable 
film-converter camera as shown. The zone plate has a series of equi- 
area annular zones which are alternately transparent and opaque to 
thermal neutrons and consist, for example, of aluminum and cadmium, 
respectively. Considering only a point scattering source of thermal 
neutrons, the display on the emulsion film (the coded image) is a 
shadow of the zone plate the shape and size of which depend monotoni- 
cally on the position of the source in the z-plane, (x,y); the position 
of the zone plate relative to the source plane, z^^-z; and the position 
of the camera relative to the zone plate, ^ beam of 

coherent light passes through the coded image, the light is diffracted 
into a focal spot which is an image of the original point scatterer. 
Considering any object to be a set of point scatterers, any specimen 
can be imaged. Each plane in the object will have a different focal 
plane when imaged, so that the three-dimensional information in the 
coded image can be viewed tomographically . 

B. Coded Image 

Neglecting the detector resolution*, assuming a linear film 
response, and noting the film output is different transparencies for 



*The detector resolution is idealized by considering the film and 
converter to be in the same plane and neglecting the film resolution. 
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Figure 1. Scheme for Pseudo-Hologram Formation 
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points on the film, the coded image can be represented by the neutron 
flux shape in the As shown in Figure 1 , (x**,y**) describes 

points in the film-converter plane at z = ; (x*,y*) describes points 

in the zone plate plane at z = z^; and (x,y,z) describes points in 
the scattering body. 

Let F (x,y,z,f^), the scattering rate density, represent the 
s — 

number of neutrons per unit time per unit volume per solid angle 
scattered in the volume dxdydz about x,y,z and in the direction d^ 
about Let G(x,y,z,^ x" ,y** ,Z2 represent the angular flux 
Green* s function to the film-converter plane, then 



-n(r,r") 



G(x,y,z,£ -> x",y",Z2,^") = 



r-r 






where and ^** are vectors from the origin (see Figure 1 ) to points 

in the scattering body (x,y,z), to points in the zone plate plane 
(x*,y*,z*), and to points in the film-converter plane (x** ,y** ,Z2)> 
respectively; and ri(^,r_**) is the **optical path length** between r_ and ^**. 

In this case it is convenient to consider ri(r^,r^”) in two parts, 
the *'optical path length** in the scattering body, 9 and the 

'*optical path length** outside the body, Then 

( 2 ) 

where the second exponential on the right hand side represents the zone 
plate transparency which is defined by a set of step functions or a 
series expansion^. 
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The neutron flux, $(x",y",Z 2 ), at the film-converter plane is 



then 



(3) 



$(x",y”,Z2)= 




If an angular emitted neutron current from the surface of the body can 



the zone plate removed, = 0 in equations (3) and (4) above, 

and the uninterrupted scattered neutron flux at (x",y") in the film- 
converter plane is described. 

To evaluate equations (3) or (4) an expression for the emitted 
neutron current or the scattered neutron source must be developed via 
a suitable neutron transport approximation. The accuracy of the source 
function can be verified by computing the uninterrupted scattered flux 
from equations (3) or (4) and comparing its shape with measured values. 
As noted above, the flux shape in the (x",y") plane is sufficient to 
represent the coded image. 



be described as K (x,y,z,J^), then 
s ““ 



(4) 





^ S 



where represents integration over the body surface, dr is a surface 



element of the body and in this case = 0. Note that with 
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III. SCATTERED FLUX FROM A REGULAR SOLID 

For the development of a coded image of the interior structure to 

be practicable, the thermal neutron scattering cross section, E , of the 

s 

scattering object must be much greater than the thermal neutron absorp- 
tion cross section, E • This suggests the possible application of the 

a 

diffusion approximation to neutron transport in developing a scattered 
source function as used in equations (3) or (4), where some approximation 
will be required to define the directional dependence of the scattered 
source function. 

A. Scattered Source Function 

It is assumed that an adequate test case is to examine the 

scattered flux from a regular solid in a beam of thermal neutrons 

of intensity (neutrons per unit time per unit area) as shown in 

Figure 2, where the composition of the solid is homogeneous and such 

that E » E . As mentioned earlier, either a scattering source function 
s a 

or an angular emitted current source function can be used in computing 
the scattered flux. To develop a scattering source function, the flux 
in the solid, $(x,y,z), is derived by diffusion theory methods, and 
the solid is considered an isotropic volumetric source of strength 
s 

4tt . The scattered flux at some point away from the solid, 

(Xo>yo>Zo) in Figure 2 would be the uncollided contribution from each 
point in the solid. An angular emitted current source function on the 
surface of the solid is defined directly from the emitted current 
derived via diffusion theory, where some emission angular dependence 
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Figure 2, Test Case Geometry Used in Determining the Scattered Flux 
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must be assumed. The use of diffusion theory in computing reflection 
coefficients and the relative ease in computing the scattered flux from 
an angular emitted current source suggests the use of the angular 
emitted source, 

The diffusion equation to be solved is 



where D is the diffusion coefficient of the solid and S(x,y,z) is an 
isotropic volumetric source. It is convenient, in terms of describing 
suitable boundary conditions, to let S(x,y,z) represent the first 
collision scattering density, i,e.. 



is assumed that the first collision scattering density is isotropic. 
The boundary conditions are that $(x,y,z) is zero at the extrapolated 
boundaries of the solid. 

The solution to equation (5) (see Appendix A) obtained by the 
method of eigenfunctions for the geometry in Figure 2 is 




(5) 




( 6 ) 



where is the total thermal neutron cross section of the solid and it 



$(x,y,z) = I sin 



Imn 



Ittx 



sin 

c 




(7) 
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where 



32 Z I K 

hrm - [^J ^ 




L = D/Z , B = b+4D, c = C+4D, a = a+4D 
a 




Noting the symmetry in Figure 2, x = c/2 and y = a/2 is the center 

of symmetry in the xy-plane, and considering points (x ,y ,z ) which do 

o o o 

not see the ends of the solid, i.e., 2D<z <b+2D, then the emitted 

o 

current need only be defined for faces 1 and 2 of the solid. The emitted 
current from face 1, (c+2D,y,z), by Fick’s Law is 



/ (c+2D,y,z) - - f 



94>(x,y,z) 

9x 



( 8 ) 



X = C+2D 



where 



a^CxjY.z) 


I is- A, cos 


sin \^' 


8x 


X = C+2D Imn ^ ^ L ^ - 


L a J 



sin 
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and the emitted current from face 2 is 



(x,a+2D,z) 



^(x,a+2D,z) 

4 



D 

2 



8$(x,y,z) 

9y 



y = a+2D 



( 9 ) 



where 



8$(x,y,z) 

8y 



= a+2D 





B. The Scattered Flux 

To use equation (4), with r|(^,r^*) = 0, in calculating the 

uninterrupted scattered flux, $(x ,y ,z ) at (x ,y ,z ), it is necessary 

o o o o o o 

to define some angular dependence to the emitted current. Realizing 
it is difficult to make any generalization about the angular dependence 
of the emitted current, two possibilities for defining an angular 
dependence in terms of isotropy are: 

1. Replace the emitting face with a plane semi-isotropic source 
(hereafter termed the semi-isotropic emitted current assumption) of 
strength j"**(x,y,z) defined by equation (8) or (9). 

2, Assume the current density, j^(x,y,z,^), is semi-isotropic 
(hereafter termed the semi-isotropic current density assumption) 
then 



J 



+ 



(x,y,z) 



II. j (x,y,z,^)d^ 

2tt 



( 10 ) 
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where j(x,y,z,^) = CQ and ^ is the unit vector normal to the emitting 
area. 



Based on the semi-isotropic emitted current assumption the angular 
emitted current from face 1, K^^(c+2D,y,z,ft) , is 



J (c+2D,y,z) 

Kgj^(c+2D,y,z,Q) = — 



( 11 ) 



and the angular emitted current from face 2, , is 



K ^(x,a+2D,z,Q) 

S2 — 



(x,a+2D,z) 

2tt 



( 12 ) 



Based on the semi-isotropic current density assumption, the 
angular emitted current from face 1, K^ 2 (^“^ 2 D,y,z,^), is 

Ksi(c+2D,y,z,Q) = n,. j (c+2D ,y ,z ,^) = cos (j) (13) 



where cos == n. Q and (|) is then the angle between the neutron direction 
and unit normal to the emitting face. Evaluation of equation (10) 
yields 



C 



1 



J^(c+2D,y,z) 

7T 



(14) 



Similarly the angular emitted current from face 2, K^2(^»^"*-2D,z ,fi) , is 



K _(x,a+2D,z,fJ) = 
sz — 



J (x,a+2D-z) 

y ^ 1 



7T 



cos (j) 



(15) 
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Using equation (4), with f|(r^,r^") = 0, the uninterrupted flux, 

^’(x ,y ,z ) at (x ,y ,z ), recalling that 2D<z <b+2D, is determined 
o o o o o o o 

for the two angular emitted currents as follows: 

1. For the semi-isotropic emitted current assumption: 

a. If 2D<y <a+2D and x >c+2D then only face 1 contributes 
^o o 

to the flux at (x^,y^,z^) and evaluation of equation (4) over 
and yields 



a+2D b+2D 






J (c+2D,y,z)dydz 

X 



2D 2D 2tt J[x^-(c+2D)]^+(y^-y)^+(z^-z)^J 



(16) 



b. Similarly, if 2D<x <c+2D and y >a4-2D then only face 2 

o o 

contributes to the flux at (x ,y ,z ) and 

o o o 



C+2D b+2D 






J^(x,a+2D , z)dydz 



2D 2D 2tt j^(x^-x) [y^"(a+2D)2 ^+(z^-z)^J 



(17) 



c. If y >a+2D and x >c+2D, then faces 1 and 2 contribute 
o o 

to the flux at (x ,y ,z ) and 
o o o 



$(x ,y ,z ) = equation (17) + equation (16) (18) 

o o o 



2. For the semi-isotropic current density assumption: 

a. If only face 1 contributes to the flux at (x ,y ,z ) the 

o o o 

evaluation of equations (4) over Q and yields 



m t 



« 




m 



-■xT.s^r- 
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a+2D b+2D 



<^>(x ,y ,z ) = 
o o o 



cos (j)^J^(c+2D,y,z)dydz 
2D 2D 7T |^jx^-(c+2D)] ^+(y^~y)^+(z^“z)^J 



(19) 



where the integration of equation (4) over Q requires that 



cos cj)^ = cos (() 



r’»-r 
r=n 



— |r”-r| 



X -(c+2D) 
o 



|^[x^-(c+2D)] ^+(y^-y)^+(z^-z)^ J 



2 -] 1/2 



b. If only face 2 contributes to the flux at (x ,y ,z ) then 

o o o 



C+2D b+2D 






cos (f)^J^(x,a+2D,z)dydz 

2D 2D 7T j(x^-x)^+[^y^-(a+2D)J ^+(z^-z)^J 



( 20 ) 



where 



cos (f)^ = 



y - (a+2D) 
o 



|^(x^-x) +|y^-(a+2D)^ +(z^-z) j 



211/2 



c. If faces 1 and 2 contribute to the flux at (x ,y ,z ) then 

o o o 



$(x ,y ,z )= equation (19) + equation (20) 
o o o 



( 21 ) 



To evaluate the above integrals, a numerical technique is 
conveniently applied. The series in equations (8) and (9) are uniformly 
convergent so that the series can be integrated term by term. Thus each 
term in equations (16), (17), (19) and (20) can be integrated numerically 



15 



and a suitable convergence established. Results of scattered fluxes 
for a specific regular solid computed as outlined above are contained 
in the next section. 
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IV. METHOD OF MEASUREMENT AND CALCULATION 
A. Method of Measurement 

Measurements of the scattered flux were obtained with an arrangement 

simulating the geometry of Figure 2. As shown in Figure 3, a regular 

solid specimen, a polyethylene (C^H^) block (2.86 cm square x 5.16 cm 

long), was placed in a beam of thermal neutrons. A BF^ detector was 

used to count scattered neutrons. The detector longitudinal axis 

remained perpendicular to face 2 and in the yz~plane at x = 1.43 cm + 

2D, where D is the diffusion coefficient in cm. The center of the 

detector window was located at the following points (x ,y ,z ) for 

o o o 

measurement : 

1. (1.43 cm + 2D, 7.86 cm H- 2D,z ) where z was incremented by 

o o 

by .5 cm from .5 cm + 2D to 4.5 cm H- 2D. 

2. Same as 1, except y^ = 12.86 cm + 2D. 

The Pennsylvania State University TRIGA Reactor (Breazeale 

Nuclear Reactor Facility) provided the source neutrons. A D 2 O tank 

between the beam port entrance and the reactor core provided sufficient 

thermalization such that the thermal flux at the beam port exit at a 

7 2 

reactor power of lOOKW was on the order of 10 neutrons per cm per 
second. Collimation was achieved by placing in the beam port two iron 
plugs with 1-1/2 inch holes followed by a copper plug with a 3/4 inch 
hole. Smaller collimations would reduce background scattering, however, 
to maintain the specimen totally within the beam the above sizes were 
the minimum allowable. The specimen was placed 15 inches from the beam 
port exit and at this point, it was estimated that the beam had diverged 
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to approximately twice the exit diameter with a flux on the order of 
6 2 

10 neutrons per cm per second. The cadmium ratio in the beam was 
approximately 50. 

The BF^ detector was 14.5 inches long with a 1 inch outer diameter 
and was filled to a pressure of 60 cm Hg with BF^ 96% enriched in B^^. 
The detector was shielded with a cadmium jacket, and the window had a 
cadmium cover with a 1/16 inch diameter hole. By placing a neutron 
shutter between the beam port exit and the detector, an average signal 
to noise ratio of 2 was achieved. As shown in Figure 3 the neutron 
counting arrangement was typical for a gas filled proportional counter. 

The polyethylene specimen was held in place by an aluminum clamp 
attached to a horizontal bar 1/4 inches in diameter. Three of the 
side faces were covered with a cadmium sheet which reduced background 
scatter from the aluminum holding stand, and maintained the integrity 
of the non-reentrant boundaries. 

One minute background counts were taken at each point described 
above with only the specimen removed. The gross counts were also 
taken for one minute and were sufficiently high to yield counting 
fractional standard deviations on the order of .025. The counts 
obtained for points 5 cm 7.86 cm + 2D) and 10 cm (y^+ 12.86 cm + 

2D) from face 2 were normalized with the z^ - 2D = 2.5 cm point. The 
results are displayed on Figures 4 and 5. As discussed earlier, the 
normalized curve represents the relative scattered flux shape as 
amended by the detector response. 

Additional measurements were made at 2.5 cm (y^“ 5.36 cm + 2D) 
from face 2. The measurements were repeated with various strips of 



cadmium on face 2 as follows: 
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Figure 4. Calculated Flux Shape vs. Measured Flux Shape at Points (1. 43cnrf2D, 12.86cnH-2D 
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Figure 5. Calculated Flux Shape vso Measured Flux Shape at Points (1. 43cnd-2D, 7.86cnrl-2D 
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1. A 1 cm wide strip beginning at x = 1.5 cm + 2D and ending at 
X = 2.5 cm + 2D. 

2. A 2 cm wide strip beginning at x = 1.5 cm + 2D and ending at 
X = 3.5 cm + 2D. 

All counts obtained were normalized with count obtained at (1.43 cm 
+ 2D, 5.36 cm + 2D, 2.5 cm + 2D). The results are displayed on Figures 
6 , 7 , and 8 . 

B. Method of Calculation 

The method of calculation was outlined in the previous section, 
and in this case equations (17) and (20) applied. To compare the 
calculations with the measured values, a response function was inserted 
representing the BF^ detector response. The response function was 
determined as follows. 

Neglecting the finite area of the detector face, the probability 
of detection for a neutron entering the detector at an angle 9 with 
the detector face is 



— B 

1 - exp(-E^ r/cos 6) 



( 22 ) 



-B -B 10 

where Z = Z is the thermal neutron absorption cross section of B 
t a 

at 60 cm Hg, r is the radius of the detector and r/cos 9 represents 
the path length of the neutron in the detector. For the geometry in 
Figure 3 



cos 9 



(x -x) ^+(z -z)^ 
o o 



(x^-x) + jy^-(a+2D^ +(z^-z) 



1/2 



(23) 
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Figure 6. Calculated Flux Shape vs. Measured Flux Shape at Points (1.43cm+2D, 5.36cnH-2D 
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Figure 7. Calculated Flux Shape vs. Measured Flux at Points (1.43cm+2D, 5.36cm+2D 
With a 1cm Cadmium Strip on Emitting Face 
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Figure 8. Calculated Flux Shape vs. Measured Flux at Points (1.43cnH“2D, 5.36cm4“2D, 
With a 2cm Cadmium Strip on Emitting Face 
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Considering the finite area of the detector window, a unit flux 
in the 9 direction to the detector face produces a current of sin 9 
entering the detector, where from equation (23) 



sin 9 



[y^-(a+2D)]^ 

(x^-x) ^+^y^-(a+2D)^ ^+(z^-z)^ 



1/2 



(24) 



The detector response function, R(x^ ,y^,z^,x,y ,z) , is then the fraction 
of the flux seen, equation (24), times the probability of detection, 
equation (22), 

Assuming a semi-isotropic emitted current and accounting for the 
detector response, equation (17) was modified to 



C+2D b+2D 



o o o 



R(x ,y ,z ,x,y,z) J (x,a+2D, z)dxdz 

0 ,0 o V ’ 

2D 2D 2tt j^(x^-x)^+|y^-(a+2Di^+(z^-z)^ J 



(25) 



Assuming a semi-isotropic current density equation (20) was modified to 



C+2D b+2D 



(26) 



$(x ,y ,z ) = 
o o o 



R(x ,y ,z ,x,y,z)cos 9^ J (x,a+2D,z)dxdz 

OOP ^ y 

2D 2D T\ |(x^-x)V[y^-(a+2D)] V(z^-z)^J 



Equations (25) and (26) were solved, as noted previously, by integrating 
each term numerically using as IBM system 370 Mod 168 computer and 
Argonne National Laboratory subroutine DROMB^^, The series was truncated 
when the higher mode integral terms were less than .1% of the existing 
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sum. The thermal neutron diffusion properties of polyethylene were 
obtained from the ANL-5800^^. Calculations were made for each point 
measured by the detector and were similarly normalized. The results 
for the semi-isotropic emitted current and for the semi-isotropic 
current density are displayed on Figures 4 and 5. 

The calculations were repeated for the detector at 2.5 cm from 
face 2 as used in obtaining the measured values. Calculations were 
made at each point measured by the detector and were similarly 
normalized. Calculations with the cadmium strips on face 2 were made 
the semi-isotropic current density assumption only. The results are 
displayed on Figures 6, 7, and 8. 
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V. ANALYSIS OF RESULTS 

In analyzing the results on Figures 4 through 8 it is helpful to 
bear in mind that: 

1. The detector had a forward counting preference which decreased 
as the detector was moved away from the emitting surface (see equations 
(22) and (24)). 

2. The emitted current, used in determining the angular emitted 
current, resulted from a diffusion calculation. 

3. The angular emitted current was determined from assumptions 
regarding angular dependence. 

As shown by Figure 4, at 10 cm from the emitting face the measured 
flux shape was generally well represented by the calculated flux shapes 
with differences only at points near the boundaries of the emitting 
face. Coupled with the detector forward counting preference, the 
boundary errors would be expected from a diffusion calculation. The 
differences, as displayed in Figure 4, indicate that the boundary 
errors were the result of calculated emitted currents which were low 
near the edges of the emitting face. 

Additionally, a portion of the observed differences can be 

attributed to the assumptions used in defining an angular emitted 

current as the calculated flux shape for the semi-isotropic current 

emission differs from the calculated flux shape for the semi-isotropic 

current density emission. To quantitate somewhat, the maximum 

difference between the calculated relative flux for a semi-isotropic 

emitted current and the measured relative flux occurred at z - 2D = 

o 

4.5 cm and was 15.1%. The difference between the calculated relative 







I 
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flux for a semi-isotropic current density and the measured relative 
flux at the same point was 18.3%. Thus, the errors displayed by 
Figure 4 were due primarily to errors in the calculated emitted current 
near the boundaries. 

As shown by Figure 5, when the detector was 5 cm from the emitting 
face, the boundary errors increased. This would be expected as the 
detector preference for points directly in front increased and the 
calculated emitted current was in error at the boundaries. Moreover, 
errors in the assumed angular dependence of the emitted current 
increased. To quantitate again at - 2D = 4.5 cm, the difference 
between the calculated relative flux for a semi-isotropic emitted 
current and the measured relative flux was 26.8%. The difference 
between the calculated relative flux for a semi-isotropic current 
density and the measured relative flux was 33.8%. 

As shown by Figure 6, when the detector was 2.5 cm from the 
emitting face, the differences increased further as would be expected. 
The differences in the calculated values for the semi-isotropic current 
density and the semi-isotropic current were small compared to the total 
errors. The principal source of the errors was again the calculated 
emitted current. Moreover, it should be noted that the block dimensions 
were only IL x 2L, where L represents the diffusion length. 

The measured and calculated flux shapes at 2.5 cm from the emitting 
face with 1 cm and 2 cm cadmium strips on the emitting face are 
displayed on Figures 7 and 8, respectively. The boundary errors dis-- 
played by Figure 6 are not as evident here, due to the normalization. 

In this case a significant portion of the face which had a reasonably 



29 



correct calculated emitted current was eliminated by the cadmium strip. 
As a result, the calculated flux at the - 2D = 2.5 cm point had a 
larger contribution from the boundary areas yielding a low value. 

When normalized to this low value, the low values near the edges of 
the emitting face were dampened, decreasing the errors in the 
calculated relative flux. 

The above comments can be summarized as follows. Differences 
between the calculated values and the measured values were due primarily 
to boundary errors in the calculated emitted current, with very 
minor effects due to the arbitrary assumptions used in defining the 
angular dependence to the emitted current. As the detector forward 
preference was decreased, the observed differences decreased and the 
calculated flux shapes represented well the measured values. Although 
large differences existed between the calculated and measured values 
when the detector had a strong forward preference, measured dips in the 
relative flux over a strong absorbing strip were represented by the 



calculated values. 
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VI. CONCLUSIONS 



A. Suinmary 

The diffusion approximation was applied in calculating emitted 
neutron currents from a homogeneous regular solid in a beam of thermal 
neutrons. The calculated emitted current, with an assumed angular 
dependence, was used in computing the scattered flux shape about the 
solid. Calculated results, modified for the response of a cylindrical 
gas filled proportional counter, were compared to values measured with 
a BF^ proportional counter at points which received flux contributions 
from only one face of the solid. The observed differences in the 
calculated and measured flux shapes were due primarily to errors in 
the calculated emitted current about the edges of the face, with very 
minor effects due to the assumed angular dependence of the emitted 
current. As the directional preference of the detector decreased, 
the observed differences decreased. Moreover, regardless of a signifi- 
cant directional detector preference, trends in flux shapes over strong 
absorbing strips were well represented. It can be concluded that for a 
detector with a lesser directional preference, such as a film-converter 
detector, the calculated flux shape would better represent the 
measured flux and would be sufficient to complete the analytical 
representation. 

® • Suggestions for Future Work 



The description of the scattered neutron flux shape represents 
the first step in an analytical representation of the pseudo-hologram. 
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Future work is suggested to complete the analytical representation by 

including the response functions for the Fresnel zone plate previously 

described and for a t}rpical film-converter detector. The detector 

response function will have to be suitably approximated as was done 

herein for the BF^ detector response function. The flux shapes so 

determined could then be compared to calculated flux shapes based on 

an ideal detector and an ideal Fresnel zone plate in an effort to 

develop a correction factor for pseudo-holograms practically obtained. 

One possible application of a correction factor would be in correcting 

12 

pseudo-holograms which are decoded numerically on a computer system 
The correction factor, depending on a general application, would at 
least decrease some of the resolution loss experienced in the real 
system. 



32 



APPENDIX 

SOLUTION OF DIFFUSION EQUATION 

The geometry is as shown in Figure 2, Then the extrapolated 
boundaries are 

c = c + 4D 
a = a + 4D 
b = b + 4D 

where c, a and b are the block dimensions. 

The equation to be solved is 



DV $(x,y,z) - Z ^>(x,y,z) + S(x,y,z) = 0 
a 



(A-l) 



where the source is 



S(x,y,z) = / 



Z I e ~ 2 j^(z- 2D) 
so for 2D<z<b-2D 



0 otherwise 



(A-2) 



and the boundary conditions are ^(x,y,z) extrapolated _ q 

boundary 

To solve (A-l) use 



(A-3) 



V^4>(x,y,z) + B^^^(f)(x,y,z) = 0 



(A-4) 



to generate a linear independent, complete, orthogonal set of functions 
{4>lmn(x,y,z) } for l,m,n = 1,2,3,**»* where 
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^Imn 



s in (-^) s in s in (-^) 

cab 



and 



(A-5) 



Imn 



(i!L) + (gl) + fe' 

cab 



(A-6) 



Then 



S(x,y,z) = I 



Imn 



(A-7) 



4(x,y,») - Z Aj__^ 



Imn 



(A-8) 



From the orthogonality property of {({) i^^ (x,y ,z) } and equations 
(A-2) and (A-7), then 



32Z I K 
s o 



Imn 2 fp2 /nir. 2l 

blmTT + (^) J 



(A-9) 



where 



K = e 



-Z.b 



• TTn (6-20) 1 _ nri^ 






fniT (6-20)1 

L 6 J 



F , /nTT20^ , niT /mr20^ 

* 1" <=os(-5— ) 



(A-IO) 



Substituting equations (A-7) and (A-8) into equation (A-1) , using 
equation (A-4) and the linear independence of {c()^^(x,y ,z)} , then 



Imn 



Z (1+L^B? ) 

a Imn 



(A-ll) 
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where 



^} = D/Z 

a 

Then using equations (A-8) , (A-9) , (A-10) and (A-11) 



$(x,y,z) 



= I 



Imn 



*lmn 



(A-12) 
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